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We consider systems whose steady-states exhibit a nonequilibrium phase transition from an active 
state to one -among an infinite number- absorbing state, as some control parameter is varied across a 
threshold value. The pair contact process, stochastic fixed-energy sandpiles, activated random walks 
and many other cellular automata or reaction-diffusion processes are covered by our analysis. We 
argue that the upper critical dimension below which anomalous fluctuation driven scaling appears 
is dc — 6, in contrast to a widespread belief (see |l| for an overview). We provide the exponents 
governing the critical behavior close to or at the transition point to first order in an e = 6 — d 
expansion. 
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Absorbing- state phase transitions. Nonequilibrium 
phase transitions between an active and an absorbing 
state are encountered in a variety of fields ranging from 
chemical kinetics to the spreading of computer viruses . 
From a theoretical standpoint absorbing state transitions 
form natural counterparts to equilibrium phase transi- 
tions. The transition rates used in the stochastic dynam- 
ics employed to model the physical phenomenon under 
consideration do not satisfy detailed balance (with re- 
spect to an a priori defined energy function). In spite 
of this apparent freedom, the number of universality 
classes that the transition can fall into is incredibly small. 
Among known universality classes, that of directed per- 
colation (DP) is by far the broadest. And indeed, in the 
absence of additional symmetries or conservation laws, 
as was conjectured twenty years ago by Grassberger 
an absorbing state transition will invariably fall into the 
DP universality class. Only if the dynamics possess addi- 
tional features such as particle number parity conserva- 
tion, or the coupling to an auxiliary field (whether con- 
served or not, static or diffusing, etc.), or the existence of 
an infinite number of absorbing states, will the transition 
belong to a new class. Apart from the formal interest in 
classifying nonequilibrium phase transitions and in iden- 
tifying which microscopic ingredients make a transition 
belong to a given universality class, there exists a greater 
challenge. As recently summarized by Hinrichsen [Q, in 
spite of its domination in the theoretical physics litera- 
ture, the directed percolation univerality class was actu- 
ally never observed in a real experiment. This is often 
attributed to the presence of defects in a real experiment 
(DP is known to be very sensitive to quenched disorder) 
or to other ill-controlled effect, such as a hidden conser- 
vation law or a coupling to an auxiliary field. 

On one hand the interest in absorbing state transi- 
tions was recently enhanced by the discovery by Vespig- 



nani, Zapperi and coworkers |5|] of their relationship with 
the ubiquitous phenomenon of self-organized criticality 
(SOC) It was estabhshed that the scahng behav- 
ior observed there was entirely governed by an under- 
lying phase transition (thereby, incidentally, tempering 
the mystics of SOC). On the other hand, in a separate 
wave of articles, research has focused on absorbing-state 
transitions in which the order parameter freezes into one 
among an infinity of absorbing states, but without any 
additional conservation law. The paradigmatic example 
of a system showing such a behavior is the pair contact 
process, initially introduced by Jensen and Dickman [Q, 
for which Muhoz and coworkers [|| devised a convincing 
phenomenological picture that we shall later use as our 
starting point. 

The existence of an infinite number of absorbing 
states (in the large-system limit) and the coupling to an 
auxiliary static field are the common characteristics to 
the microscopic models that we wish to investigate here. 
It is the purpose of the present work to provide a full 
renormalization group picture of the phase transition 
at work in systems possessing an infinite number of 
absorbing states, with or without an additional con- 
servation law. We shall rely on a phenomenological 
Langevin approach as a starting point. From there 
we shall show how renormalization group arguments 
can be applied. This will lead us to finding the upper 
critical dimension of those models. Then we will sketch 
the reasoning leading to the expression of the critical 
exponents within the framework of an expansion around 
the upper critical dimension. As will be clear there is 
interesting new physics to learn from the many techni- 
cal obstacles that pave the way to the full scaling picture. 

We now turn to a presentation of two microscopic 
models chosen for their representativeness and ease 
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of formulation in which the absorbing state phase 
transitions we wish to study appear. We also introduce 
the Langevin equations encoding their dynamics. Then 
we sketch the field-theoretic line of reasoning leading to 
the computation of the critical exponents. Finally we 
provide a critical discussion of our results in the light of 
the existing literature. 

Two models and their phenomenological coarse-grained 
description. In the Pair Contact Process (PCP), parti- 
cles are thrown on a lattice. A pair may either annihi- 
late A -I- ^4 ^ or produce a single offspring A + A ^ 
A + A + A. Each lattice site can be occupied by at most 
one particle. Since isolated particles cannot diffuse, a 
configuration in which pairs of nearest neighbors are ab- 
sent remains frozen in time. The order parameter -0 is 
local density of pairs of nearest neighbor particles and 
the control parameter is the branching over annihilation 
rate ratio. The Langevin equation believed to describe 
the dynamics of ip was coined by Mufioz et al [|| . It reads 

PCP: dti^iv, t) ^DAMr, t) - a^{v, t) - .9iV''(r, t) 
-g3V'(r,i)e--fod*V(r,t')+^(j.^i) 

(1) 

with T] a Gaussian white noise whose correlations are 
(?7(r,i)77(r',i')) = g2#^''Hr " " i')- The coefH- 

cients ct, gi , 32 and 53 are coarse-grained analogs of the 
reaction rates. The memory term is the signature of 
the feedback of isolated particles on the pair dynamics. 
For a detailed explanation of the physical origin of the 
various contributions appearing in Eq. (|l]) we refer the 
reader to The list of microscopic models whose 

coarse-grained description is believed to be encoded in 
the Langevin equation (|l|) also includes the Transfer 
Threshold Process ll and various models for catalysis 



(the dimer reaction 
dimer-trimer models 



the dimer-dimer or the 
|). A detailed analysis of the 



mean-field properties of this equation was provided by 
Munoz et al. 

The second family of models that we are interested in 
is embodied by the so-called Manna or stochastic fixed- 
energy sandpiles (FES) |l^, Grains of energy are 
initially randomly distributed on a lattice. Whenever a 
lattice site is occupied by more than 2 particles (in di- 
mension d) the excess particles randomly hop to a nearest 
neighbor site. The number of active sites plays the role 
of the order parameter. The total number of particles is 
strictly conserved. An evolution equation for the local 
density of active sites ij) was recently proposed by Dick- 
man, Vespignani et al. It reads 

FES: dti^ = Dl\,^{v, t) ~ a^{v, t) - gi^b^{r, t) 

/•* (2) 
+ gi4>{r,t) At' A,i;{r,t')+ri{v,t) 
Jo 



with ?7's correlations having the same expression as in 
the PCP case and the coefficients 17,(71,(72,54 are coarse- 
grained analogs of the microscopic transition rates which 
depend on the conserved quantity. The nonlocal memory 
term expresses that space fluctuations of the static and 
conserved field have a feed-back effect upon the order 
parameter dynamics. Again a rich variety of models 
were shown to be described by the same coarse-grained 
dynamics, such as the Conserved Transfer Threshold 
Process, Activated Random Walkers [Q, or some model 
for epidemic spreading in which healthy individuals are 
static 

There are some hazards in relying on naively built 
phenomenological equations as many other interactions 
are generated by a coarse-graining procedure. In princi- 
ple, all terms allowed by symmetries should be included. 
We shall later see that, indeed, some relevant symmetry- 
preserving terms have to be considered. Let us now re- 
call the common mean-field behavior of those models. 
Denoting for both models by a the deviation of the con- 
trol parameter from its critical threshold value, in the 
steady-state the order parameter behaves as 



0" 



= 0, CT > 



(3) 



In mean- field we simply have ^ — —ajgx. At the critical 
point, the order parameter decays according to 



il)(t) ^ t^\ t large 



(4) 



The correlation length and the relaxation time diverge 
as (T ^ according to f ^ |c| " and r ~ \a\^'^^, 
respectively. Within the framework of a mean-field 
analysis one finds the following values for the critical 
exponents: (3 — 1,6 = 1,1^ = ^ and z = 2. 



Analytic strategy. The interplay between short-time and 
short-distance fluctuations with long range correlations 
lies at the heart of the anomalous scaling observed in 
the vicinity of a critical point ("anomalous" to be un- 
derstood as non mean-field). Short distance singularities 
(usually referred to as ultraviolet divergencies) govern 
the way scaling properties deviate from their mean-field 
expressions. The usual strategy is to retain the leading 
UV divergencies and to perfrom a renormalization 
group analysis with those only. In the vast majority 
of cases, this is sufficient to reach physical conclusions. 
Sometimes however, contributions that were thrown 
away in the course of the analysis are crucial in preserv- 
ing the correct physics, though they are irrelevant in 
determining the renormalization group fixed point. This 
is the present situation. 

Let us start with the PCP. The picture is the following. 
We expand the exponential memory term and find that 
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the contribution g^tp J* ^ exhibits the leading short-time 
and short-distance singularity. Power counting shows 
that the bare dimension of 3253 is d — 6 (in units of a 
length) . This type of interaction is known to describe the 
dynamical percolation universality class |l^, which 
has an upper critical dimension dc = 6. Yet, in the 
present case, the phase diagram would not be repro- 
duced correctly if only the truncated expansion of the 
exponential were kept. As seen on the mean-field evolu- 
tion equation, and as already noted many times in the 
literature the asymptotics are governed by the local 
nonlinear term —gi^"^ in Eq. (|^). However 51(72 has bare 
dimension d — 4 which signals that it is in fact a dan- 
gerously irrelevant coupling, along with the subsequent 
powers of the argument of the exponential. The cou- 
pling gi is not relevant in determining the underlying 
fixed point structure but essential in preserving the over- 
all phase diagram. An important consequence is that the 
usual scaling assumption for the order parameter 

{4;{t)) =b-'^Tib-H,b^/''\a\), Marge (5) 

must be abandoned since the scaling function J- will ex- 
hibit a singular behavior in the gi variable as the latter 
approaches 0. In mean-field J- depends on gi as l/gi- 
While irrelevant variables are traditionnally omitted from 
the list of arguments of J-, in the present case this would 
lead to unphysical conclusions. The next step to see how 
this dangerously irrelevant coupling gi is renormalized to 
in the vicinity of the dynamical percolation fixed point 
(a necessity overlooked in [^). In order to do this we have 
followed the technical procedure recalled by Janssen and 
Schmittmann (see also the references therein). We 
skip all technical steps and merely quote the final 
result: gi{b) ~ as the coarse-graining scale b is 

increased. We have denoted by £ = 6 — d and the result 
for the exponent is given to first order in e. In order 
to obtain the leading s behavior we have combined the 
mean-field expression for the scaling function J-' appear- 
ing m Eq. (|) as far as its gi{b) dependence is concerned 
with the scaling properties of the field, time and a at the 
dynamical percolation fixed-point. This has led us to the 
following critical exponents, 

/3 = l-^£, -^-l-ie, (6) 

the expressions of which are given to leading order in 
e = 6 — d. 

As far as the FES described by Eq. (^ are concerned, 
the situation is a bit more delicate. A rather involved 
RG analysis ||T^ shows that in fact the dynamical 
percolation vertices are generated already at one-loop 
order (the one-loop graph obtained by connecting 32 
and two g^s leads to effective gi and 33 vertices, and 
no effective 34). This is because, as is very often the 
case with gradient interactions, the short-scale behavior 



of the memory term —gi%l;{v,t) &t' ^r'4'{v,t') converts 

into an effective —gzik{v,t) j^dt'il^^Tji') contribution 
(and higher powers) after coarse-graining. This also 
explains the failure of naive power counting directly 
on (H) which, naively performed, would lead to the 
erroneous conclusion that dc = 4 since 3254 has bare 
dimension d — 4 while 3253 has dimension d — 6. And 
then we may apply the reasoning of the previous para- 
graph and the results of Eq. (||) continue to hold for FES. 

What the field-theory approach teaches us can be 
summarized as follows. First, a detailed analysis of the 
renormalized interactions shows that both classes of 
models considered are described by the same-field theory 
which has an upper-critical dimension dc = 6. This 
result is in contradiction with the existing literature of 
the last ten years . Critical exponents are expected to 
behave differently from mean-field in space dimensions 
d < 6. We have determined the critical exponents of 
the phase transition to leading order in e = 6 — d. As 
to FES-like systems, at the renormalization flow fixed 
point, the conservation law is irrelevant. This seemingly 
innocuous property leads to technical difficulties since 
the conservation law is crucial in saving the overall phase 
diagram. In both the PCP and the FES cases, while the 
underlying fixed point is that of dynamical percolation, 
computing the critical exponents requires to perform an 
independent renormalization of an irrelevant coupling. 
Thus the critical exponents cannot be deduced from 
those of percolation by means of hyperscaling relations. 
The validity of our approach holds in the vicinity of 
dc — Q. As d is lowered, already at d = 4 new qualitative 
features show up in the theory. Since the percolation 
fixed point becomes trivial in d = 1 an entirely new 
picture must inevitably set in as dimension is lowered, 
perhaps as early as d is decreased below d = 4 and 
possibly again at d = 2. Splitting of the universality 
classes of the PCP and FES in low space dimensions 
must be envisaged as well. 

There are several criticisms that can be opposed to the 
present work. First of all, the use of phenomenological 
Langevin equations sometimes proves hazardous, but, 
fortunately other formalisms |^ allow exact mappings 
onto field-theories having exactly the same features as 
the ones discussed in this work. Indeed, the effective 
Langevin equation (||) used for sandpiles is, strictly 
speaking, not correct, since coarse-graining (loop correc- 
tions) shows that other dominant terms are generated. 
Most importantly, however, we must admit that there 
is absolutely no numerical evidence supporting our 
findings. An obviously too concise summary of the 
numerical state of affairs is that both FES and 
PCP belong to the DP class in dimension one (recent 
studies disagree with the results of |2^, such as p^). 
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But as dimension is increased PCP is still found to 
belong to the DP class while FES are found to form an 
independent universality class. Liibeck, in some recent 
simulations claims that dc = 4 by performing 

simulations directly in high dimensions. A cheap way 
out for the theoretician is to refer the reader to a recent 
preprint by Grassberger [ p3| in the somewhat different 
but related context of forest-fires in which dynamical 
percolation plays some role as well. There numerical 
proof is provided that it might well be impossible with 
present day computers to ever reach the true asymptotic 
scaling regime. We recall that there are many obstacles 
on the numerical side: the impossibility to use simple 
finite-size scaling relations (due to the presence of a 
dangerously irrelevant operator). Quenched disorder 
(for FES |25|| ) and other long-term memories are known 
to be difficult to overcome numerically, and finally the 
proximity of the directed percolation fixed point whose 
infiuence must be felt until the system eventually crosses 
over to its actual asymptotics. All of those features, 
we fear, play a part in rendering the reaching of the 
true asymptotics a hopeless endeavor. As far as the 
observed lack of universality of spreading exponents is 
concerned we believe that the phenomenon can be 
understood within a renormalization group picture. As 
demonstrated in a much simpler case in |27 spreading 
exponents exhibit nonuniversal values that depend on 
the moments of the initial distribution of both the order 
parameter and auxiliary fields. This dependence is all 
the stronger as the initial distribution deviates from a 
strict Poissonian law. In less favorable cases the initial 
distribution introduces couplings at the initial time 
that cannot even be renormalized, thus questioning the 
existence of a scaling regime. 

As a conclusion we summarize our findings. We 
have provided field-theoretic arguments showing that a 
large number of stochastic processes exhibiting a phase 
transition between an active state and an absorbing 
state (picked up among an infinite number of such) 
belong to the universality class of dynamical percolation 
and therefore have = 6 as their upper critical dimen- 
sion (instead of dc = 4 as appears in the last decade 
literature). The existence of a static field to which the 
order parameter is coupled is the common feature to all 
considered models. We have provided expressions for the 
critical exponents within the framework of renormalized 
perturbation theory in the vicinity of the upper critical 
dimension to leading order in e ~ Q ~ d. A technical 
account of the results presented in this letter is in 
preparation [ p^ . 

The author akcknowledges several discussions with HK 
Janssen, HJ Hilhorst, A Vespignani, R. Pastor-Satorras, 
M.A. Mufioz, H. Chate and J. Kockelkoren. 



[1 
[2: 

[s: 

[4] 
[5] 

[6: 

[7 

[s: 

[9 

[10: 
[11 

[12 

[13 

[14 

[15 

[16 
[17 

[18 

[19 

[20 



[21 
[22' 
[23^ 
[24' 



[25; 



[26' 
[27; 



Frederic. van- Wi jland® th.u-psud.fi 

R. Dickman, cond-mat 0110043, Invited talk to Statphys 
XXIII. 

J. Marro and R. Dickman, Nonequilibnum phase tran- 
sitions in lattice models (Cambridge University Press, 
Cambridge, 1999) ; H. Hinrichsen, Nonequilibnum Crit- 
ical Phenomena and Phase Transitions into Absorbing 
States, Advances in Physics 49, 815 (2000). 
P. Grassberger and A. de la Torre, Ann. Phys. (N.Y.) 
122, 373 (1979). 

H. Hinrichsen, Braz. J. Phys. 30, 69 (2000). 
R. Dickman, A. Vespignani, and S. Zapperi, Phys. Rev. 
E 57, 5095 (1998); A. Vespignani, R. Dickman, M. A. 
Mufioz, and S. Zapperi, Phys. Rev. Lett. 81, 5676 (1998). 

H. J. Jensen, Self-Orgamzed Criticality (Cambridge Uni- 
versity Press, Cambridge, 1998). 

I. Jensen and R. Dickman, Phys. Rev. E 48, 1710 (1993). 
M. A. Mufioz, G. Grinstein, R. Dickman and R. Livi, 
Phys. Rev. Lett. 76, 451 (1996) and Physica D 103, 485 
(1997). 

M. A. Mufioz, G. Grinstein and R. Dickman, J. Stat. 
Phys. 91, 541 (1998). 

J. F. F. Mendes, R. Dickman, M. Henkel, and M. C. 

Marques, J. Phys. A: Math. Gen., 27 3019 (1994). 

E. V. Albano, J. Phys. A: Math. Gen. 25, 2557 (1992); 

I. Jensen, Phys. Rev. Lett. 70, 1465 (1993). 

J. Kohler and D. ben-Avraham, J. Phys. A 24, L621 

(1991). 

A. Vespignani, R. Dickman, M. A. Mufioz, and S. Zap - 



peri, Phys. Rev. E 62, 4564 (2000); |gond-mat/0003285 . 
S. S. Manna, J. Phys. A 24, L363 (1991); D. Dhar, Phys- 
ica A 263, 4 (1999). 

M. Rossi, R. Pastor-Satorras, and A. Vespignani, Phys. 
Rev. Lett. 85, 1803 (2000). 

R. Pastor-Satorras and A. Vespignani, Phys. Rev. E 62, 
R5875 (2000) and Eur. Phys. J. B 19, 583 (2001). 
H.K. Janssen, Z. Phys. B 58 311 (1985). 
J.L. Cardy, J. Phys. A 16, L709 (1983) and P. Grass- 
berger and J.L. Cardy, J. Phys. A 18, 267 (1985). 
H.K. Janssen and B. Schmittmann, Z. Phys. B 64 503 
(1986). 

R. Pastor-Satorras, A. Vespignani and F. van Wijland, 
in preparation. 

M. Doi, J. Phys. A9, 1465 (1976); L. Peliti, J. Physique 
46, 1469 (1985); B. P. Lee and J. Cardy, J. Stat. Phys. 
A 80, 971 (1995). 

S. Liibeck, Phys. Rev. E 64, 016123 (2001). 
J. Kockelkoren and H. C hate, in preparation . 



P. Grassberger, preprint cond-mat /0202022[ . 
R. Dickman,M. Alava, M.A. Muoz, J. Peltola, A. Vespig- 
nani and S. Zapperi, Phys. Rev. E 6 4, 056104 (2001); R . 
Dickman, T. Tome, M.J. de Oliveira, [sond-mat/0203565. 



The initial static field fluctuations play the same role 
as quenched disorder added to a simple DP process [ p6[ , 
which cannot easily be accounted for in a Langevin equa- 
tion for -i/) only. 

H.K. Janssen, Phys. Rev. E 55 6253 (1997). 

F. van Wijland, K. Oerding, and H. J. Hilhorst, Physica 

A 251, 179 (1998). 



